Degeneracy between mass and spin in black-hole-binary waveforms 



Emily Baird, 1 Stephen Fairhurst, 1 Mark Hannam, 1 and Patricia Murphy 1 

School of Physics and Astronomy, Cardiff University, Cardiff, CF24 3AA, United Kingdom 

(Dated: November 6, 2012) 

We explore the degeneracy between mass and spin in gravitational waveforms emitted by black-hole binary 
coalescences. We focus on spin-aligned waveforms and obtain our results using phenomenological models 
that were tuned to numerical-relativity simulations. A degeneracy is known for low-mass binaries (particu- 
larly neutron-star binaries), where gravitational-wave detectors are sensitive to only the inspiral phase, and the 
waveform can be modelled by post-Newtonian theory. Here, we consider black-hole binaries, where detectors 
will also be sensitive to the merger and ringdown, and demonstrate that the degeneracy persists across a broad 
mass range. At low masses, the degeneracy is between mass ratio and total spin, with chirp mass accurately 
determined. At higher masses, the degeneracy persists but is not so clearly characterised by constant chirp mass 
as the merger and ringdown become more significant. We consider the importance of this degeneracy both for 
performing searches (including searches where only non-spinning templates are used) and in parameter extrac- 
tion from observed systems. We compare observational capabilities between the early (~2015) and final (2018 
onwards) versions of the Advanced LIGO detector. 



I. INTRODUCTION 

The advanced LIGO and Virgo detectors are likely to allow 
us to observe numerous black-hole-binary coalescences in the 
coming years [T] |2) . While the detectors are being installed, 
the challenge is to devise search and source extraction meth- 
ods that will identify all black-hole-binary signals in the data 
whilst also allowing for the most accurate estimation of the 
source parameters possible. 

Our understanding of gravitational waveforms emitted by 
black-hole mergers has improved dramatcially over the past 
five years, with the aid of a large number of numerical simula- 
tions, from which several models of the waveforms have been 
constructed (see e.g., Ref [3| for an overview). A large frac- 
tion of the simulations have focussed on non-precessing sys- 
tems (where the black holes are either non-spinning or have 
spins aligned with the orbital angular momentum); see, for 
example, Refs. H and 0. Other than for high black-hole 
spins (greater than a ~ 0.8), this space has been rather well 
covered for comparable-mass binaries, and several waveform 
models are available i^HHl. The models characterize wave- 
forms based on the masses of the two components and, in the 
case of the models we will consider in this paper, the total spin 
of the system. We will make use of these waveform mod- 
els to investigate degeneracies in the waveforms. Since we 
are only considering a three-dimensional subspace of the full 
eight-dimensional space of binary masses and spins, many ef- 
fects, most notably precession, cannot be probed. However, 
it is interesting to consider this smaller space because it is 
very likely that any degeneracies will persist in the full pa- 
rameter space. Furthermore, recent results have shown that 
the waveform for a precessing binary can be factorized into 
precessional effects and a non-precessional part that is well 
modelled by the three parameters we consider here . 

In post-Newtonian (PN) theory, there is a well known de- 
generacy between the black holes' mass ratio and the total 
spin, which arises at 1.5PN order, while a combination of the 
total mass and mass ratio (the "chirp mass") remains relatively 
well constrained ifTolfTTl . Here, our focus is on higher-mass 
waveforms that include a merger and ringdown. We demon- 



strate that the degeneracy persists in the full waveforms and 
also up to higher masses where only the later part of the inspi- 
ral and the merger/ringdown are in the detector's sensitivity 
band. However, for the high-mass binaries, the degeneracy 
in mass and spin is not so clearly characterized by the chirp 
mass. 

A degeneracy in the emitted gravitational waveform across 
the parameter space has numerous consequences for gravita- 
tional wave (GW) searches, both good and bad. The posi- 
tive effect is that a degeneracy in the parameter space will 
reduce the volume that needs to be searched, thereby reduc- 
ing the computational cost and trials factor associated with the 
search. The majority of searches of LIGO and Virgo data have 
made use of non-spinning components in the template wave- 
forms (see e.g., Refs. lfl2l[T3ll ). A degeneracy between mass 
and spin means that the template waveforms have covered a 
larger fraction of the parameter space than might be naively 
expected. Furthermore, it has recently been argued that a two- 
dimensional bank of templates is sufficient to cover the space 
of spinning neutron-star binaries [ 14 1. In this paper, we inves- 
tigate the effect of using non-spinning waveforms in a search 
for black-hole binaries with spins and show that the mass-spin 
degeneracy renders the search more sensitive to spinning sys- 
tems than might be expected. Nevertheless, a search using 
waveforms incorporating spin effects would be a significant 
improvement and the degeneracy should make it computation- 
ally feasible. 

On the other hand, a degeneracy has a negative effect on 
parameter estimation. To make the most of gravitational- 
wave observations, accurate extraction of the physical param- 
eters is of paramount importance. There are detailed multi- 
dimensional methods under development to accurately re- 
cover the parameters l[T5l4l8l . However, there's nothing that 
can be done about a real degeneracy in the emitted waveforms 
— there is no way of telling them apart. We evaluate the ef- 
fects of mass-spin waveform degeneracy on our ability to ac- 
curately recover masses and spins and discuss the astrophysi- 
cal implications. In the process we introduce a simple method 
to estimate the parameter-estimation confidence intervals. 

Throughout the paper, we will provide sample results us- 
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ing a number of waveforms at different masses and mass ra- 
tios. We also make use of a number of noise curves for the 
advanced detectors, to illustrate how this effect is likely to 
change as the detectors approach their final, design sensitiv- 
ity. No attempt has been made to perform an exhaustive study 
across the full parameter space, and this is left as a future 
project. 

The layout of the paper is as follows: in Sec. [11] we outline 
the waveform models that we use, and discuss some of the 
assumptions that went into producing them and their range of 
applicability. We also discuss the interpretation of results in 
terms of mismatches between waveforms and introduce the 
detector noise curves used in our studies. In Sec. [Ill] we focus 
on the degeneracy between mass ratio and spin for low-mass 
binaries, giving the post-Newtonian argument for this degen- 
eracy and illustrating the degeneracy for a number of cases 
and different detector sensitivities, and discuss implications 
for searches. Then, in Sec.[lV]we extend the results to higher 
mass binaries and show that the degeneracy persists, although 
in a different form. In Sec.[V]we discuss implications for the 
accurate estimation of parameters and astrophysical inference. 

II. MODEL AND METHODS 

A. The phenomenological waveform models 

We describe the GW signal from black-hole binaries with 
non-precessing spins (i.e., the spins are aligned or anti-aligned 
with the binary's orbital angular momentum) using the phe- 
nomenological models presented in Refs. [6| and [7|. For 
consistency with the labelling used within the LIGO-Virgo 
Collaboration |[T9l we refer to these models respectively as 
"PhenomB" and "PhenomC". ("PhenomA" refers to an earlier 
model of non-spinning binaries IT2014221 .) In both models the 
waveforms are parametrized by their total mass M = m\ +1112, 
mass ratio 77 = m\ni2/M 2 , and an effective total spin parame- 
ter, 

1 / x 
X = ~^\ m ^X\ +m 2 X2) 

= Xs + SXa, (1) 

where %i = Sj/mj for each black hole with angular momentum 
S(, 8 = (mi —m2)/M, and the symmetric and anti-symmetric 
combinations of the spins are % s = (%\ +X2)/2 an d Xa = 
(Xi ~ /fc)/2. The phenomenological models incorporate a 
PN description of the inspiral, while the merger and ringdown 
regimes are tuned using the results of numerical simulations. 

Both models have the same basic structure. The waveform 
is represented in the Fourier domain as h(f) = A(/)e !>1 '^. 
The amplitude A (/) and phase are modelled separately, 
using input from PN theory (inspiral), the observed properties 
of NR waveforms (plunge-merger), or results from perturba- 
tion theory (ringdown). The models are all power series in 
the frequency /, and the coefficients in the model are written 
as polynomials in the two physical parameters 77 and % (the 
total mass is an overall scale factor), and it is the coefficients 
of these polynomials that are then calibrated to hybrids of PN 



and NR waveforms. For both models, the amplitude is con- 
structed in a similar manner, with expressions for each of the 
inspiral, plunge-merger, and ringdown portions of the wave- 
form modelled independently. In PhenomB O the three parts 
are connected as piecewise functions, and in PhenomC |7| 
smooth tanh-function interpolation is used. To obtain an ex- 
pression for the phasing, PhenomB uses a single series expan- 
sion, matching the coefficients beyond leading order to hy- 
brid waveforms, as well as results from the test mass (tj — > 0) 
limit. PhenomC, on the other hand, uses the complete Tay- 
lorF2 PN inspiral phasing, and only the late inspiral/merger 
phase is fitted in a narrow frequency range to numerical sim- 
ulations, while the ringdown waveform is obtained from an- 
alytically derived quasi-normal mode expressions for the fre- 
quency and attached continuously to the merger phase. There 
are 54 free parameters in PhenomB, and 45 in PhenomC, al- 
though the final models are both functions of only {M, r),x}- 
There are three notable differences between the two models: 
(1) the PhenomB PN-NR hybrids are produced in the time do- 
main, using TaylorTl for the PN part, and the PhenomC hy- 
brids are produced in the frequency domain, using TaylorF2 
for the PN part [23 1; (2) PhenomB incorporates information 
about the test-mass limit; (3) in PhenomC the phase evolution 
during inspiral incorporates PN calculations up to 3.5PN or- 
der (although the spin terms are complete only up to 2.5PN), 
while in PhenomB only the leading-order PN inspiral term is 
fixed, and the remaining terms up to 3.5PN order are tuned to 
the hybrid waveforms. 

There is good qualitative agreement between the two mod- 
els 0, although no detailed quantitative comparison has yet 
been performed. We choose to work with PhenomB for gen- 
erating the results presented in this paper. We have also cross- 
checked some of the calculations against PhenomC, and we 
comment further on this in Sec. lIVI 



B. Detectors and noise curves 

In this paper we compare spinning- and non-spinning- 
binary signals with reference to the expected sensitivity of the 
Advanced LIGO detector (aLIGO) I24jj26l . The sensitivity 
curves we use are shown in Fig. [T] During the early science 
runs, expected around 2015, the advanced LIGO detector is 
unlikely to be at its full design sensitivity. Consequently, we 
use the "early aLIGO" noise curve lP27ll to give results indica- 
tive of what may be achieved in the early runs. At its opti- 
mum sensitivity several years later, the anticipated sensitivity 
is given by the "zero-detuned high-power" noise curve 11281 . 
In this paper we will take that as the "final" design sensitivity 
of the detector. Over the parameter space of binaries that we 
study in this paper, the early aLIGO curve represents a sen- 
sitivity of roughly five times greater in the signal-to-noise ra- 
tio (SNR) than in the initial LIGO detectors during their final 
S6 science run, which corresponds to an increase in potential 
sources of two orders of magnitude. The final curve repre- 
sents a further factor of three improvement in SNR, or about 
30 times as many potential sources over early aLIGO. Finally, 
we also consider the "no signal recycling" (no-SRM) config- 
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1 — x 2 /2. Thus, the match is directly related to 
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FIG. 1 . Noise curves for the advanced LIGO detector configurations 
that we consider in this paper: "early", "no-SRM" and "final", which 
corresponds to the zero-detuned, high-power configuration. See text 
for more details. 



uration of the detector, that could be achieved by the aLIGO 
detector operating without a signal recycling cavity. This has 
comparable low-frequency sensitivity to the final configura- 
tion, but significantly worse sensitivity at high frequencies. 
Although it is unlikely to be an observational mode, the non- 
signal-recycled curve offers a means to compare the effects of 
low and high frequency sensitivity upon our results. 



C. Waveform Mismatches 

We use the standard inner product between two waveforms, 
h\(f) and %(/) with respect to the power spectral density 
S n (/) of a detector | 



(h\h 2 )=4Re 



' *!(/)*£(/) 
Sn(f) 



df. 



(2) 



The match between two waveforms is defined as their normal- 
ized inner product, maximized over time and phase shifts of 
the waveform, 



M(/.a 2 ) = max^M. 

A(,A0 |«l||«2| 



(3) 



Alternatively, we can consider the situation of a true wave- 
form produced by a physical source, hj, and a model wave- 
form Iim that we will use to search for the real signal in de- 
tector data. If we normalize both waveforms, then we can 
split the model waveform into a component parallel to the 
true waveform, plus a component that is orthogonal to the true 
waveform, in the sense of our inner product. In other words, 
we have 



the (relative) amplitude of the "error" waveform. 

In a GW search, a template bank of model waveforms is 
constructed [ 30 1 such that the match between every point in 
the waveform parameter space and the nearest template in the 
bank is at least 0.97. Assuming that the model waveforms are 
physically correct, such a template bank ensures that we will 
lose no more than 10% of signals in our search. (A match of 
0.97 means that the sensitivity range of the detector is only 
97% of its optimum, and the detector is therefore sensitive to 
only 0.97 3 ps 0.9 of its optimum volume, and so we lose about 
10% of signals.) 

If the physical waveforms do not agree exactly with the 
model waveforms, this will lead to an additional loss in 
match between a signal and the best matched template, and 
consequently a reduction in the number of signals observed 
above threshold. In addition, to counter non-stationary detec- 
tor noise, a number of signal consistency tests are included 
into analyses to distinguish signals from non-stationarities 
or "glitches" in the data Oil . These tests are used to ei- 
ther remove completely any transients that do not match the 
templates or else to down-weight their significance. Since 
searches performed to date have made use of non-spinning 
waveform families these thresholds have been set relatively 
loosely (and tested with spinning signals) to ensure they are 
not removing signals. For the high matches between signal 
and template we consider in this paper, the effect of signal 
consistency tests will be minimal, and we will not consider it 
further. 

In the next sections we will identify the regions of the non- 
spinning waveform parameter space that provide a match of 
greater than 0.97 with the chosen waveform (which will usu- 
ally incorporate spin). In doing so, we densely sample the 
non-spinning waveforms to identify all points with a match 
above 0.97. When performing a search, there will then be 
two contributions to the mismatch between the signal and the 
best matched template: one due to the difference between the 
waveform and search space and a second arising from the dis- 
crete sampling of the template space. The match between the 
signal and the closest template is guaranteed to be above 0.94 
as the mismatches add linearly in this case (see e.g. ||32l for 
details). We are therefore requiring that the potential loss of 
SNR due to a mismatch between the model waveform and true 
waveform is no greater than the maximum possible loss due 
to the discreteness of the template bank. 

We will also investigate how the parameters of the best- 
match non-spinning waveform vary when the signal corre- 
sponds to a non-precessing binary. In Sec. |V]we will see how 
the match can be used to give an estimate of the parameter 
estimation accuracy. 



III. DEGENERACY BETWEEN r] AND % AT LOW MASSES 



\J 1 —x 2 hj +xh[. 



(4) 



where hp is the normalized "error" waveform that satis- 
fies (/irl^fi) = 0. In this form, the match is MQitJim) — 



We first consider the degeneracy between the symmetric 
mass ratio 77 and the effective total spin % of the binary. This 
effect is already well known in PN theory ifTOiriTl . which we 
will discuss first, and then we will look at inspiral-merger- 



ringdown (IMR) signals at three different values of the total 
mass, M = {20,50, 100}M o . 



A. Degeneracy in PN theory 

The phase evolution of a compact binary in PN theory has 
been calculated up to 3.5PN order in the non-spinning terms, 
and up to 2.5PN in the spin effects; see Ref. [23 1 for a recent 
summary of PN treatments of the phase. Up to the leading 
order that includes spin, the phase for non-precessing binaries 
is given in the frequency domain by 
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where v = (ftM/) 1 / 3 . If we define the chirp mass of the bi- 
nary as ^£ = Mr] 3 / 5 , then we see that the leading factor is 
proportional to 7f/) 5 / 3 , and the phase evolution is dom- 

inated by the chirp mass. This motivates the observation that 
in GW searches we expect to measure the chirp mass with 



high accuracy, and we will see examples of this in Sec. Ill B 



At the next-to-leading order, the phase evolution depends on 
the mass ratio 77, and dependence on the spins enters at the fol- 
lowing order. We can absorb all of the spin effects at this order 
into an effective spin term, Jpn = Xs + ^Xa ~ (7677 / 1 13)x?, 
and this is what is proposed in the inspiral template family 
discussed in Ref. [33 1 (and this is proportional to the leading- 
order spin-orbit parameter j3 used in Refs. ifTOlfTTI ). The phe- 
nomenological models (6) UJ use the simpler effective spin 
% = Xs + &Xa to describe full inspiral, merger and ringdown 
waveforms. 

If we adopt for the moment the PN effective spin term, we 
can write the deviation from the leading-order term as 
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Thus, it is possible to mimic the effect of spin by modify- 
ing the mass ratio (while keeping the chirp mass constant). 
It is clear, however, that the required modification to 77 will 
vary as v changes over the course of the chirp, as expected as 
this is only an approximate degeneracy. The majority of the 
power from a binary merger is accumulated between around 
30 to 300 Hz, and over this range v changes by only a factor 
of two. This explains why the approximation is reasonable 
across the whole signal, and this leads to the common claim 
that non-spinning templates can be used to detect GW signals 
from spinning binaries, but there will be an offset in the mea- 
surement of the mass ratio and the total mass. It may also 
be possible to exploit this degeneracy to search for spinning 
binaries using a non-spinning model, but with 77 extended to 
unphysical values, 77 > 0.25, in order to cover more of the 
spinning-binary parameter space; a similar idea has already 
been suggested to extend inspiral-only searches beyond their 
expected region of validity 11341 . It is the degeneracy between 
mass ratio and spin that we will investigate here. 
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FIG. 2. Match between spinning 2OM waveforms with non- 
spinning templates in early aLIGO. Each of the strips shows the re- 
gion of M-77 space for which the non-spinning waveform has a match 
of 0.97 or higher with a spinning signal. The total spin is indicated 
next to each region. For spins above 0.5 or below -0.7 the best over- 
lap with a non-spinning waveform is less than 0.97. The curve of 
constant chirp mass is indicated by a dotted line. 



B. Degeneracy in full IMR models 

We would now like to investigate whether the degeneracy 
from post-Newtonian theory is present in full IMR waveform 
models. 

We consider a set of spinning signals and those non- 
spinning waveforms (i.e. the phenomenological model with 
X — 0) that provide a good match to the signal. The efficacy 
of a waveform model in a GW search can be estimated by cal- 
culating the best match (fitting factor) between any member 
of the waveform model, and the target signal. Therefore, this 
will give an immediate indication of the merits (or otherwise) 
of using a non-spinning waveform model to search for merg- 
ers of spinning black holes. We consider fitting factors above 
0.97 to be adequate as discussed previously. 

Fig. [2] shows the results for a binary with a total mass of 
2OM and mass ratio 1:4 (TJ = 0.16), i.e. a4M -16M binary, 
and matches calculated using the early aLIGO spectrum. We 
consider sources with a total spin in the range [-1,1] in steps of 
0.1. For each configuration, we identify non-spinning wave- 
forms which give a match of greater than 0.97 with the spin- 
ning source. When the binary spin is zero, the non-spinning 
model matches the target signal at the correct parameters (in- 
dicated by a star), and along a strip in the parameter space 
with a width of 5% in mass, and 10% in 77. When the bi- 
nary contains spinning black holes, the non-spinning model 
matches the signal, but with a bias in the mass and mass ratio. 
The binary spins are indicated with different colours, labeled 
by the total spin of the binary. For large-spin signals there 
are no non-spinning waveforms that have a match above 0.97. 
When using the early aLIGO noise curve, the range of spins 
for which non-spinning waveforms have a match above 0.97 
is [—0.7,0.5]. However, for spins above about X = 0.3, the 
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FIG. 3. Same results as in Fig. [2] but now parametrized in terms 
of chirp mass and symmetric mass ratio. This figure reinforces the 
accuracy with which the chirp mass can be recovered. 



best-match waveform has an unphysical value of the mass ra- 
tio tj. The curve of constant chirp mass Mtj 3 / 5 is indicated by 
a dotted line. 

The first thing to note is that, in all cases, the chirp mass 
of the best matched non-spinning waveform is very similar to 
the true chirp mass of the system. We can see from Fig. |2]that 
the high-match regions follow a line of constant chirp mass. 
In Fig.[3]we re-parameterize the results of Fig.|2]in terms of 77 
and ^ . We see that, even though there is a strong degeneracy 
between 7} and %, the correct chirp mass ^# = 6.66M Q falls 
within the high-match region for most values of spin. The 
deviation from the correct chirp mass is no more than ~ 3% 
for large anti-aligned spins. 

The second important point relates to the mass-spin degen- 
eracy. For aligned-spin binaries (% > 0), the best match is 
obtained with a lower-mass, higher-mass-ratio model, while 
for anti-aligned spins, the best matched non-spinning wave- 
form has a higher mass and lower mass-ratio. This is to be 
expected as aligned spins are expected to cause the system to 
"hang up" and chirp more slowly, mimicking a lower-mass 
signal. 

It is informative to return to the post-Newtonian phasing 
given in Eq. |5]) to see whether the observed degeneracy 
matches what is theoretically expected. In low-mass binaries, 
the signal is dominated by the inspiral, which can be repre- 



sented in a PN expansion, as described in Sec. Ill A The chirp 
mass is determined to high accuracy by the leading-order term 
in the PN expansion of the phase in the Fourier domain, and 
we can then solve Eq. |5]l (with fixed chirp mass) to find the 
symmetric mass ratio that mimics the effect of the spin. We 
do this by solving |5]) as a function of frequency (in the detec- 
tor's sensitive band) and then averaging over the values of rj 
that we obtain. Fig.|4]shows the symmetric mass ratio that cor- 
responds to each value of the spin, as predicted from Eq. d3), 
and as found in the mismatch analysis of the phenomenologi- 
cal models that is shown in Fig. [2] We see that the PN estimate 
is remarkably close to that found for the full IMR models, at 




FIG. 4. Comparison of the %-r\ degeneracy as found for the phe- 
nomenological model (stars), and from the requirement that = 
A < ¥(xpn = 0) in Eq. (|5} for fixed chirp mass, for signals with 
M = 20Mq and 77 = 0.16. The naive PN prediction agrees well with 
the full IMR results for low spins, and for most anti-aligned spins. 



least for the 20 M e binaries used in this example. 



C. Evolution of noise curves 

We present results for two design aLIGO configurations in 
Fig. [5] the no-signal-recycling (no-SRM) configuration and 
the final high-power zero-detuned sensitivity curve (final). 
Qualitatively, the results are quite similar for the no-SRM and 
final noise curves. Again spinning signals are recovered with 
a match above 0.97 using non-spinning templates, with a de- 
crease in recovered mass for positive spin systems. The chirp 
mass is still well recovered. However, the detector's extra sen- 
sitivity makes it easier to distinguish the spinning signals with 
a non-spinning model. The range of spins for which we can 
find a non-spinning signal with matches greater than 0.97 is 
now only % E [-0.4,0.3] for no-SRM and % € [-0.3,0.3] for 
the final configuration. We also see that each match region 
shrinks, although its location is unchanged. In moving from 
the early noise curve to no-SRM, the most significant sensi- 
tivity improvement is at low frequency, while the final con- 
figuration offers much greater high frequency sensitivity. The 
results for these two curves are quite comparable (and sig- 
nificantly better than the early curve) indicating that it is low 
frequency sensitivity that provides the biggest improvement 
in breaking the degeneracy. With the later noise curves, the 
variation of chirp mass between the spinning signal and non- 
spinning template is less than 1%. Q 



1 The figures may seem to indicate that a 21) 17 . binary with total spin of 0.5 
that would be detected by a non-spinning search in the early aLIGO con- 
figuration, but not later ones. This is of course not true: the final detector 
is roughly three times more sensitive at these masses, and so a signal with 
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FIG. 5. Matches of spinning 2OM waveforms with non-spinning 
templates in aLIGO at design sensitivity. The upper panel shows the 
results for the no-SRM configuration, while the bottom shows results 
for the final sensitivity 



IV. MASS-SPIN DEGENERACY AT HIGHER MASSES 

For black-hole binaries with masses greater than 2OM the 
merger and ringdown parts of the signal become increasingly 
important and contribute an ever increasing fraction of the 
signal-to-noise ratio. At higher masses, we do not expect the 
chirp mass to determine the waveform to such an extent as for 
the 20M© system and there is no a priori reason to expect that 
a degeneracy between mass and spin will persist. 

We begin by considering a 50M e binary with mass ratio 
1:4. At this mass, the merger and ringdown will provide a 
significant fraction of the signal to noise ratio. As a crude es- 
timate of the effect, imagine that the inspiral part of the wave- 
form is valid up to the innermost stable circular orbit (ISCO) 
of a point particle orbiting a Schwarzschild black hole, which 
is at 90 Hz for a 50M Q binary. For the early aLIGO noise 
curve, the inclusion of the merger and ringdown will increase 



SNR 10 in early aLIGO would have an SNR of 30 in the final configura- 
tion. Even with a match of 0.9, this would still give an SNR of 27, and the 
match would be sufficient that it would pass any signal consistency tests. 
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FIG. 6. Matches of spinning 50Mq waveforms of mass ratio 4:1 
with non-spinning templates in early and final aLIGO. Each of the 
strips shows the region of M-rj space for which the non-spinning 
waveform has a match of 0.97 or higher with a spinning signal. 



the SNR by a factor of 4. At the final aLIGO sensitivity, the 
improved low frequency response increases the contribution 
of the inspiral but the merger and ringdown still contribute a 
comparable SNR to the inspiral (so that the overall SNR is 
~ 50% higher, because they add in quadrature). Thus, there is 
no reason to expect that the chirp mass is still well recovered 
or that the degeneracy discussed previously still holds. 

In Fig. [6] we again show the regions of the non-spinning 
parameter space which give a match greater than 0.97 with 
a spinning binary (with values of % from -1 to 1 in steps of 
0.1). Here, we show the results for both the early and final 
aLIGO spectra. The results are qualitatively similar to the 
2OM binary with spinning signals well recovered by the non- 
spinning waveforms, and the chirp mass accurately recovered. 
The sizes of the regions are roughly consistent with the lower 
mass system, with a mass accuracy of ~ 10% and range in 
tj of 0.05-0.1 for the early noise curve. For both early and 
design curves, signals with spins between 0.2 and -0.4 have 
matches above 0.97 with non-spinning waveforms. Interest- 
ingly the degeneracy still roughly follows the line of constant 
chirp mass, even though the merger and ringdown contributed 
significantly to the SNR of the signal. 

Next, we increase the mass to 100M e and repeat the anal- 
ysis. At this mass, the point-particle ISCO is at 40 Hz, so 
there is essentially no power from the inspiral in the initial 
detectors, and only a small amount in the early aLIGO. Even 
at the final sensitivity, the merger and ringdown provide the 
vast majority of the SNR (a factor of 4 more than the inspi- 
ral). Thus one expects that it is really the merger that dictates 
the waveform as seen by the detector. Figure [7J shows results 
for 1OOM , again with r\ = 0.16. For aLIGO at final design 
sensitivity, the curves are again lying pretty much along the 
line of constant chirp mass. 

We do not show results for the early aLIGO sensitivity 
curve in Fig. [7J because the results depend strongly on the 
choice of model, either PhenomB or PhenomC. The variation 
of the two models with respect to the physical parameters is 
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Binary Mass (M Q ) 

FIG. 7. Matches of spinning IOOMq waveforms of mass ratio 4:1 
with non-spinning templates in design aLIGO. 



sufficiently large through merger and ringdown that they lead 
to qualitatively different results in our mismatch studies, when 
we use the early aLIGO noise curve, where the merger and 
ringdown contribute essentially all of the SNR. The two mod- 
els were developed as pioneering models of aligned-spin IMR 
waveforms for use in searches, and for this we expect that they 
are sufficient; but their fidelity with respect to parameter esti- 
mation at high masses has not yet been tested, and should be 
a focus of future work to refine them. 



V. PARAMETER RECOVERY 

Once a black hole merger has been detected, we wish to ex- 
tract the signal parameters as accurately as possible. Param- 
eter estimation proceeds by identifying regions of parameter 
space which give a signal that is most consistent with the data. 
It stands to reason that these regions will contain waveforms 
that have a good match with the observed signal. Thus we 
might expect that confidence regions in parameter estimation 
are associated with regions of high match between signal and 
template. In this section, we show this expectation to hold in 
detail in the high-SNR limit and derive an expression for the 
value of the match that corresponds to a given confidence at a 
given SNR. 

We begin by presenting the argument, using the Fisher 
matrix formulation, and then we use this connection to re- 
interpret our earlier results in terms of parameter recovery. We 
follow the formalism used Ref. [35|, but provide only a brief 
discussion of the Fisher matrix formalism, and refer the reader 
to the Ref. [ 35 1 (and references therein) for more details. 



A. Connection between mismatch and confidence regions 

Let us assume that a signal ho is present in the data. The 
detector data is given as 

s{t)=h (t)+n(t). (6) 



To investigate signal recovery and parameter extraction at 
leading order, we Taylor expand the signal in a region of the 
true parameters (0 = 0) as 

h(9) = h + 9% + ... (7) 

Here hi — d,h is used to denote the derivative of the waveform 
with respect to the paramters 0,-, which will include M, rj, % 
as well as the amplitude, phase and coalescence time that are 
maximized over in the match calculations. 
The likelihood, for a given set of 0,-, is 

, , x f (s-h(0)\s-h(0))} 

p(5|0)oc eX p|-^ UZj. (8) 

Substituting the expressions for s (|6| and h(6) (j7| into the 
above and keeping leading order terms gives 

, , , f ("I") , , ^ OiOAhAhj)) 

pfrlflJoce^j-Lli + ftfolfr)- A 2 " J> \ (9) 

In the context of Bayesian parameter estimation, this can 
be recast in terms of a posterior probability distribution for 
the parameters 0,- using Bayes theorem: 

p{0\s)<xp{s\6)p{6) (10) 

where p(9) is the prior probability distribution for the param- 
eters. In what follows, we use a uniform prior on the parame- 
ters. In general, such a prior will not be physically motivated 
but, for a detectable signal, the likelihood will be peaked in a 
small enough region of parameter space to make this approx- 
imation reasonable. 

Given the above, we are interested in calculating the ex- 
pected offset between the true parameters and the mean value 
from the posterior distribution. We also want to evaluate the 
size of a confidence region containing a given fraction p of 
the posterior probability. These quantities give us two differ- 
ent measures of the expected accuracy of parameter recovery. 

We begin by calculating the mean of the parameter 0, as 

(9i) = J dee i p{d\s) = {h i \hj)-\n\h j ). (li) 

Thus, the mean of the posterior distribution will be offset from 
the true parameter values due to the presence of noise. One 
way to characterize this is the expected size of the error wave- 
form, h£ = h((0j)) — ho as 

(h 2 E ) n = ((n\h i )(h i \h J )(hj\n)) n =k, (12) 

where k is the dimension of the parameter space and ()„ in- 
dicates the expectation value over many noise realizations. 
Thus, on average the difference between the true signal and 
"best fit" waveform will have an amplitude of \fk. 

Next, we turn our attention to confidence regions in param- 
eter space — a region © of parameter space that contains a 
given probability p of the posterior distribution, 

p=(d0p{6\s). (13) 
J® 
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There are many ways to construct such a region, and one typ- 
ically also requires the smallest possible region. To calculate 
confidence regions in the Fisher approximation, we begin by 
observing that the covariance between parameters is given as 



(didj) = (hi\hj) 



-l 



(14) 



Using this expression and Eq. ( [TTj ), we can re-express the pos- 
terior distribution as 



pi s , oc exp \ - (6,))(h,\hj)(6j - (9j)) 



: exp 



; (\h(9)-h((8))\ 7 



(15) 



Then, the minimum volume region which contains a fraction 
p of the posterior probability is the one for which 



\h(o)-h((e))\ 2 <xi(i-p), 



(16) 



where — p) is the chi-square value for which there is 
1 — p probability of obtaining that value or larger, and k de- 
notes the degrees of freedom, determined by the number pa- 
rameters included in 0,-. At leading order, the confidence in- 
terval contains all points for which the amplitude of the dif- 
ference between the model and best fit waveforms lies below 
the given threshold. 

There are six parameters in the aligned-spin waveform 
model: M, tj, %, A, Q c and t c . When calculating matches, we 
have maximized over the latter three parameters (amplitude, 
phase and time) and reported match over the three dimen- 
sional space of mass, mass ratio and spin. Thus, we have cal- 
culated the three dimensional matches M(hi\fi2), maximized 
over A, <j> c and t c . It is straightforward to re-cast our earlier 
results in terms of mismatches. To do this, we note that 



min \hy — h%\ 

A 2 



1*1 



1 



(hi\h 2 ) 2 

M 2 N 2 



(17) 



When we restrict to the subspace of masses and spins, we can 
write the final term as the match between the waveforms. 



By substituting Eq. ( 17 1 into Eq. (To I, we see that the con- 
fidence region is defined by all points in the parameter space 
for which the match satisfies 



M(h(d),h((0)) > i 



xHi-p) 

2p 2 



(18) 



where the value of k is given by the dimension of the remain- 
ing parameter space. This expression gives us a straightfor- 
ward way to re-interpret the match calculations presented ear- 
lier. For example, with a three dimensional parameter space, 
the 90% confidence region at a given SNR is given by: 



M(h(0),h((9))>l- 3 -^ 

r 



(19) 



which has the nice property that a match of 0.97 corresponds 
to a 90% confidence region at an SNR of 10. For a two di- 
mensional parameter space, the right hand side is 1—2.3/ p 2 , 



meaning that a match of 0.97 corresponds to 90% confidence 
region at an SNR of 9. 

We should note that there are numerous assumptions used 
in the derivation of these results. Most significantly, all of the 
Fisher matrix results consider only leading order effects and 
become less reliable at lower SNR, see e.g., Ref [35 1 for a 
detailed discussion of the issues. When actually calculating 
confidence regions for a signal, detailed parameter estimation 



analyses calculate the posterior distribution ( 15 i and integrate 
to find the region containing 90% of the probability. Here, we 
are using a hybrid approach: we use the result of the Fisher 
matrix calculation to decide the threshold on match required 
to define the given confidence region, but then calculate the 
match between waveforms exactly, without recourse to any 
approximations. Furthermore, we are maximizing the match 
over three dimensions, analytically for A and (j) and using a 
Fourier transform to search over time. Thus, we are only ap- 
plying the Fisher result to the three dimensional subspace of 
mass, mass ratio and spin. Even there, we are merely using 
the calculation to determine the appropriate match threshold: 
the matches are calculated exactly. Consequently, the match 
regions should be in good agreement with the 90% confidence 
regions. We verified that for an SNR of 20, the region identi- 
fied by our match criteria did contain 90% of the probability 
to within ±0.5%. 



B. Implications for detectability and parameter estimation 

Let us now return to the results of Sec. [Ill] and [IV] and 
use the relationship derived above to re-interpret the results 
in terms of confidence intervals. We can interpret Fig. [2] in 
two ways: in the context of either a non-spinning or spin- 
aligned templated search. If we were to perform a search with 
non-spinning templates and observe a M = 20M e , rj =0.16 
binary with non-spinning components with SNR 9, then the 
90% confidence region would lie along a strip in the param- 
eter space with a width of 5% in mass, and 10% in rj. For 
a binary containing the same mass black holes with spins of 
-0.5, the 90% confidence region would be of approximately 
the same size (5% in mass, and 10% in 77) but centred at 
M = 27M , 77 = 0.100 The same holds for other values of 
the spins: the reported statistical uncertainty in parameters is 
relatively small, while the systematic errors can be significant. 

Alternatively, we can consider the three dimensional space 
of M-r\-%. In that case, a match above 0.97 corresponds to 
a 90% confidence region at SNR of 10. Thus, for example, 
Fig. [2] shows that the 90% confidence region for a signal with 
m\ = 16, mi = 4 and % S [—0.7, 0.5] would contain waveforms 
with non-spinning components. Figures [3] to [7] can be inter- 
preted in a similar manner. For higher SNRs, the confidence 



2 Strictly speaking, if the best match between the signal and a model wave- 
form is less than unity, the relation between match and confidence region 
in equation ITS) must be modified to reflect this. The 90% confidence re- 
gion will contain all points with a match of 0.97 or greater with the best fit 
model waveform. 
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FIG. 8. 90% confidence intervals for a 2OM , 1:4 (77 = 0.16) sig- 
nal with spin % = 0.2, using the design aLIGO noise curve. The 
top panel shows the full three-dimensional confidence region for 
SNR=10; the gray surface indicates the % = plane, and corresponds 
to the x = 0.2 region in Fig. [5] The lower two panels show the same 
data projected onto the M-% and m\-ni2 planes, and also indicate the 
90% confidence regions for SNRs 20 (blue) and 30 (green). The true 
physical parameters are indicated by a ball (top panel), or star (lower 
panels). Regions with unphysical rj (> 0.25) are not shown. 



intervals shrink. At SNR 10 they correspond to matches of 
0.97, SNR 20 to 0.992 and SNR 30 to 0.9965. 

We can use the analysis of the previous section to estimate 
the 90% confidence intervals if we employ an aligned-spin 
model for parameter estimation. For signals with an SNR 
of 10, the three-dimensional confidence region in (M,rj,x) 
will correspond to matches greater than 0.97. In these cases 
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FIG. 9. Same as Fig. [5] but now the signal has spin % = —0.2. 



we calculate the matches between a given aligned-spin sig- 
nal, and all other aligned-spin model waveforms with varying 
mass, mass-ratio and total spin. In Fig. [8] we show the 0.97 
match volume for a 20 M & 1:4 binary with % = 0.2, using the 
final aLIGO noise curve. The top panel shows the full three- 
dimensional confidence region, and in the lower panels the 
results are projections onto the M-% and m\-m-i planes, to aid 
the interpretation. Figure [9] shows similar results, but for a 
binary with % = —0.2. 

Since our waveform model now includes spin, the best 
match is unity at the correct parameters. But from the fig- 
ures we see that the 90% confidence region extends well be- 
yond the correct parameters, and is far from the naive image 
conjured by the term "error ellipse". We note in particular 



10 



that this volume includes a region of the non-spinning binary 
parameter space; the intersection with the % — plane is con- 
sistent with the lower panel of Fig [5] If we were to estimate 
the parameters of this signal with an aligned-spin waveform 
model, we would find that the best-fit parameters indicated 
that the black holes were spinning, but could not rule out that 
they might in fact be non-spinning — or indeed have spins 
with the opposite orientation. 

The reason for this large uncertainty in the parameters is 
that the SNR is low; it is only 10. For an SNR of 20, the 90% 
confidence interval corresponds to a match of above 0.992, 
and for an SNR of 30, the match must be above 0.9965. We 
see that for these higher SNRs, the confidence region shrinks. 
However, only at SNR 30 is a non-spinning signal excluded 
from the 90% confidence region for the % = 0.2 case, and 
even at this SNR (which may be quite rare in aLIGO obser- 
vations), the smaller mass is determined to within only 25%. 
When x = —0.2 (Fig.|9| the 90% confidence region includes 
spins as high as % = —0.7, even when the SNR is 20. Note 
also that in these figures we have included only those portions 
of the confidence regions with r\ < 0.25, i.e., physically ac- 
ceptable values; the confidence region would extend further 
in the upper two panels if the model were not constrained. 



VI. CONCLUSIONS 

We have investigated the degeneracy between mass-ratio 
and spin in gravitational waveforms, going beyond inspi- 
ral models to include merger and ringdown signals from 
black-hole-binary mergers. We have used phenomenological 
inspiral-merger-ringdown (IMR) models to study the subset 
of the full binary parameter space that includes non-spinning 
black holes, and black holes with spins parallel or anti-parallel 
to the orbital angular momentum of the binary. 

A degeneracy between mass-ratio and spin is already 
known for inspiral signals, which are relevant to ground-based 
gravitational -wave detectors for masses M < 12 M |[T4l[36l . 
We find that this degeneracy persists at higher masses, where 
the detectors are also sensitive to the merger and ringdown. 
This means that in a GW search that uses only non-spinning 
binary templates (and this is computationally cheaper than in- 
cluding spin), the signal may still be detected, but the best- 
match template will have strongly biassed parameters. The 
mass-ratio-spin degeneracy follows lines of constant chirp 
mass, so the chirp mass will be recovered with reasonable 
accuracy, even up to high masses and for spinning signals. 
However, the total mass and mass-ratio of the best-matched 
template will be biased. As shown in Fig. |2j binaries with 
higher aligned spins will be recovered with a higher mass ra- 
tio, while those with anti-aligned spins will be recovered with 



a lower mass ratio. If we restrict the search to physical values 
of the symmetric mass ratio rj, then comparable-mass bina- 
ries whose spins are aligned with the orbital angular momen- 
tum will tend to be missed in the search. We show how these 
results will change as the detector evolves towards its final 
sensitivity in Fig. [5] and for higher masses in Figs. [6] and [7] 

We also demonstrated that it is possible to use match calcu- 
lations to estimate the confidence intervals in parameter esti- 
mation. For example, the 90% confidence region for the three- 
dimensional parameter space of (M, rj, %) for signals with 
SNR 10 is given by the region with matches above ^0.97. 
This allows us to estimate how the mass-rato-spin degeneracy 
will be reflected in parameter estimation. We show that for 
modest SNRs (~ 10) it may be difficult to determine whether 
a binary includes spin, and even for high SNRs (~ 30), the 
mass-ratio-spin degeneracy impairs the accurate recovery in 
the individual black-hole masses; see Figs. [8] and [9] These 
results will affect the astrophysical conclusions that can be 
drawn from GW observations, and this will be explored in 
more detail in future work. It would also be interesting to ex- 
ploit our knowledge of the mass-ratio-spin degeneracy in a 
jump proposal for Bayesian parameter-estimation codes ifTBT - 

El. 

Our results are restricted to aligned-spin binaries, and will 
be affected by the inclusion of precession effects. At present, 
no complete inspiral-merger-ringdown model for precessing 
systems exists. It seems unlikely, however, that the inclu- 
sion of precession effects will reduce the size of the confi- 
dence regions. In general it is far more likely that increas- 
ing the dimensionality of the model parameter space will in- 
crease the parameter uncertainty, because the confidence re- 
gions will now correspond to lower matches. The recent re- 
sults in Ref. [9 1 show that precession effects have only a weak 
impact on the phasing, which suggests that while the inclusion 
of precession is unlikely to break the degeneracy discussed 
here, it is also unlikely to introduce an additional degeneracy 
in these three parameters. On the other hand, the inclusion of 
higher harmonics, which were ignored in this study, may im- 
prove the accuracy of the parameters. The net impact of these 
two effects remains to be studied in more detail in future work. 
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